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The density dependent relativistic hadron field theory(DDRH) is extended to inclusions 
■ of pentaquark + .We investigate In-medium properties of + and nucleons.As well as the 

neutron and proton, the effective + mass decreases as the baryon density increases, and 
remains larger than nucleons'.The effective mass of + is calculated, Mq + ~ 0.71M@+ at 
normal nuclear density with the fraction of + 0.2. The effective masses of all baryon 
increase with the increasing fraction of + .The binding energy is also studied. We find the 
binding energy is much larger than it in nuclear matter without + ,and then discuss the 
stability of the system with conclusions of pentaquark + . It is shown that the system 
becomes bounder added small fraction of + .finally we discuss the condition when the 
system is the most stable and calculate the minimum of binding energy and the fractions of 
all baryons as functions of baryon density. 

PACS: 14.20.-c; 24.10.Nz 
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1 Introduction 

The possibility of the existence of pentaquark + (uudds)ha.s been theoretically discussed for 
many years. In 1997, Diakonov, Petrov, and Polyakov ,in the framework of the chiral soliton model, 
predicted a exotic baryon 9 + with mass M ~ 1.53Gev,S = +l,spin 1/2 and isospin jTJ . How- 
ever, the subject has been thrust to the fore- front during the past two years by the experimental 
discovery of an exotic baryon[2J El HI El El HI El HH CHI HI]. It has K + n quantum numbers(B=+l, 
Q=+l, S=+l), and its minimal quark content should be uudds. The remarkable features of 
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the + are its small mass (1540MeV) and very narrow width (<25MeV)|8. i .The discovery of the 
exotic baryon + with positive strangeness opens new possibilities of forming exotic + hyper- 
nuclei as same as in the case of negative strangeness A, E, H hypernuclei. D. Cabrera et al have 
suggested that + could be bound in nuclei and calculated self energy of the + pentaquark in 
nuclei[l2j.they shown that the in-medium renormalization of the pion in the two meson cloud of 
the + leads to a sizable attraction, enough to produce a large number of bound and narrow + 
states in nuclei. 

On other hand, the relativistic mean field approximation have been a widely used and successful 
approach to describe properties of nuclear matter and finite nuclei (T^j. the Density Dependent Rel- 
ativistic Hadron (DDRH) theory was introduced previously as an effective field theory for isospin 
nuclei fT^l I15j.ln DDRH theory the meson-baryon interactions is described by meson-nucleon ver- 
tices which are functionals of the fermion field operators. The coupling coefficients are not constants 
but field operators. DDRH theory reduces to a Hartree description with density dependent cou- 
pling constants similar to the initial proposal of Brockmann and Toki in mean-field approximation 
and its development has been achieved a significant progress in recent years [16 .For including the 
scalar isovector 5 meson, DDRH is successfully applied to asymmetric matter at extreme neutron- 
to-proton ratios [T7|. 

In the present paper we attempt to extend DDRH theory to inclusions of pentaquark + and 
investigate properties of matter with inclusions of + .We study medium modifications within the 
mean field theory and get the mass modifications of baryons In-medium. We calculate the binding 
energy at different baryon fractions. We find the system becomes more bound with inclusions of 
+ but when too much + added it becomes less bound. The critical fraction of + is 0.56 in 
isospin symmetric matter and it becomes lower in isospin symmetric matter . The isospin effects 
on the effective mass and the binding energy are also discussed .Finally, we calculate the minimum 
of binding energy and the fractions of all baryons as functions of baryon density when the system 
is the most stable. 

2 The density dependent relativistic hadron field 

The density dependent relativistic hadron field theory has been proven to give a good descrip- 
tion of nuclear matter in bulk and of the properties of finite nucle^Zj .Assuming the interaction 
between nucleons and mesons is similar to + and meson, We start from the effective lagrangian 
density which includes pentaquark + 
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i=N,e+ 

-* N r"(p)^T a * N A^ (4) 



where 



F$ = d lk A<l>-d v Af. (5) 

Here,£^ and Cm are the free baryonic and the free mesonic Lagrangians, respectively, and inter- 
actions are defined by £j ni .The meson fields are denoted by 4> a ,(fis, A w , and A p , and their masses 
by m a , m$,mu and m p , respectively. The isospin Pauli matrices are written as r a , r 3 being the 
third component of r a . For exotic baryon G + , as we assume isospin / = 0, r a = 0, there is 
no coupling to p and 5 mesons. The main difference to DDRH theory in Ref.^Tj is that the la- 
grangian density includes the interactions of pentaquark G + and mesons. In contrast to standard 
QHD models [TBI EH ED], the meson-baryon vertices f a (a = cr,uj,5,p) are not constant numbers 
but the baryon field operators \I/ dependent. 

In the mean field approximation(MFA),the meson field operators can be replaced by their 
expectation values, which are classical fields [T%1 IT§1 12U] 

4>i -> (4>i) = 4>m (6) 
A* - (A*) = Ai (7) 

The field equations are reduced to 

(-v 2 + ml)^ = rZ(p)p% + rf (p)pf, (8) 

(-V 2 +ro|)*» = I?(p)p&, (9) 

(-v 2 + ml)4 u) = ^(p)p n + rf (p)p 0+ (10) 

7 2 , 2\ Ap) _ r N{ \ 3 



(- V 2 + mjW = T>)p 3 (11) 
where the densities are values of the following ground state expectation 

Ps = (12) 

P* = (13) 

p% = ^nt-^m) =p P s -p N s (14) 

p 3 = (¥+T 8 ¥tf) = P P - P N (15) 

The Dirac equation, separated in isospin, is the only remaining operator field equation 

[i lfl d» - T^°A» - T °r 3 Ag - (Mjy - - t^*,)]** = (16) 

[i 7 „0" - rf - (M 0+ - rf^)]*o+ = (17) 
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The in-medium mass of nucleons and pentaquark + is then given by 

M P = M N -V^(p)-Vf(p)^ s (18) 
M* N = M N -T^p) + Tf(p)^ s (19) 

The field equations can be further simplified assuming translational invariance and neglecting the 
electromagnetic field. Solutions of the stationary Dirac equation 

[ 7 ^-m,K(*) = (20) 



The usual plane wave Dirac spinors is given bv[2*T] 
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where Xb is a two-component Pauli spinor and the index b distinguishes between neutrons 
, protons and pentaquark + .Due to the inclusion of the 5 meson, The effective mass differs for 
neutrons and protons. The kinetic 4-momenta /c^ M and the energy E£ of the particle are related 
by the in-medium on-shell condition, so it is easily to get k^ 2 = m^ 2 . Integrating over all states 



k < kp b inside the Fermi sphere and introducing Ep b 
in infinite nuclear matter lead to 



k Fb + ml the scalar and vector densities 
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The density-momentum tensor is defined 
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Here the meson-baryon vertices f a (a = a, u, 5, p) are dependent on the baryon field operators 
\]/ rather than constant numbers. Relativistic covariance requires that the vertices are functions 
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T a (p) of Lorenz-scalar bilinear forms p(^f, \&) of the field operators. In mean-field approximation 
they reduce to density dependent coupling. Hofmann got a rational approximation as follows by 
mapping DB calculations jTTj 
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Table 1: The DDRMF model parameter sets are taken from Ref.jTTj 



meson a 


a 




5 


P 


m a [MeV] 


550 


783 


983 


770 


a a 


13.1334 


15.1640 


19.1023 


12.8373 


b a 


0.4258 


0.3474 


1.3653 


2.4822 


c a 


0.6578 


0.5152 


2.3054 


5.8681 


d a 


0.7914 


0.5989 


0.0693 


0.3671 




0.7914 


0.5989 


0.5388 


0.3598 



p = 0.16 [fm- 3 ] 



The theta-meson coupling constants can be determined with the quark meson model developed 
in Ref. |22j.For isospin / = 0,we only consider a and p meson. The equation of motion for meson 
field operators follows 

d^a + mla = glqq, (27) 
d^tf + mlfif = g«qj u q, (28) 

where g q a and are the quark-meson coupling constants for o and u, respectively. The mean fields 
are defined as the expectation values 

(A\a(t,r)\A) = a(r), (29) 

(A\u v (t,r)\A) = 6(v,0)u(r), (30) 

Where \A) is the ground state of the nucleus and (t,r) are the coordinates in the rest frame of the 
nucleus. In the mean field approximation the sources are the sum of the sources created by each 
nucleon - the latter interacting with the meson fields. Thus 

w(t,r) = (99(^ r ))i> ( 31 ) 

i=l,A 

qYq(t,r)= £ (qYq(t,r)) u (32) 

i=l,A 

where (..■)% denotes the matrix element in the nucleon i located at Ri at time t. Finally for the 
strange baryon,X.H. Zhong etc have got the relations |23j 

„S _ n l n N v n S _ n q n N n S _ n N foo\ 
9a ~ S / B > 9 " ~ If 9 " P ~ P ' 



Where n q is the total number of valence u and d quarks in the baryon S. 

For + ,n q =4.T s / B « 1 for all hyperons in practice [21]. In our calculation We use the relations as 
follows 



rf(p) = ^( P ) (34) 



The baryon density is given by 



P = p n + P p + P 



e+ 



~ ($Kf K * Fn + (3rcy K ^ + (3tt)^ 0+ (35) 

Where Kp n ,Kp p and Kp®+ are the Fermi momenta for n,p and G + , respectively. The baryon 
fraction is 

Yi= nM P — e+ ,* = n,p,9 + (36) 

pn _|_ p P _|_ pU^ 

The energy per baryon for multi-Theta matter is defined 

§ = - - ^M n - F P M P - Y e+ M e+ (37) 
P 

3 Results and discussion 

We first calculate the effective baryonic masses with different G + fractions in isospin symmetric 
matter and isospin asymmetric matter. The results are presented in Fig. 1(a) and Fig. 1(b). It is 
shown that the effective masses all decrease as the baryon density increases,whether the nucleons or 
the B + .The isospin effects on the effective nucleon mass can not be neglected due to the interaction 
between the scalar isovector <5,especially at high density.The neutron-proton(n/p) effective mass 
splits. But the effects on the m* &+ is negligible. The effective Theta masses always decrease more 
slowly than the nucleons. The fraction Y" + can also affect on baryonic masses. At the same baryon 
densitythe baryonic masses all increase with more conclusion of pentaquark @ + . 

Fig. 2(a) shows the binding energy in isospin symmetry matter. We plot the binding energy with 
different G + fractions. Our calculation show that the binding energy is sensitive to the variation of 
the + fraction. The absolute minimum of binding energy increases until the + fraction reaches 
about 0.23,witch is up to nearly 24Mev at the baryon density 1.4po-But when the fraction of + 
reaches to 0.56, the system becomes not bound anymore. 

The binding energy of the system without proton is presented in Fig. 2(b). The system with 
inclusion of a small fraction of pentaquark + becomes stable. The absolute minimum reaches 
the peak at the fraction 0.3. More B + added, system becomes less bound. When the percent of 
pentaquark B + reaches to half an whole, system becomes unstable again. It shows that the binding 
is also sensitive to the isospin by comparing with the curves in Fig. 2(a). 

Minimizing the E/B with respect to Yq+ and Y p ,we got the minimum of the binding energy. The 
results are given in in Fig. 3. It is shown that The minimum of binding energy occurs at the density 
about 1.5po- in low-density region(0 < p < 1.5po)>the binding energy increases with the increasing 



baryon density. The opposite dependence occurs in hight-density region(p > 1.5p )- Fig. 3(b) 
displays the fractions of neutron, proton and pentaquark @ + as calculated by minimizing the 
energy. The curves display the neutron-to-proton ratio is near 1:1. It is easily to comprehend since 
the asymmetric nuclear matter is more stable than symmetric matter. 

In summarywe calculate the Multi-Theta matter could be bound in the framework of den- 
sity dependent relativistic hadron field theory.lt seems possible that there is a large fraction of 
pentaquark B + . Whether nucleons or pentaquark G + ,it's effective mass in-medium is decreasing 
with the increasing baryon density.But effective mass of pentaquark + is larger than that of the 
N, + . 
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Figure 1: effective baryonic masses. 
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Figure 3: FIG. 3(a) shows the minimum of the binding energy.FIG.3(b) shows the fractions of 
neutrons, protons and pentaquark + 
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